Entanglement control in hybrid optomechanical systems 
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We demonstrate the control of entanglement in a hybrid optomechanical system comprising an optical cav- 
ity with a mechanical end-mirror and an intracavity Bose-Einstein condensate (BEC). Pulsed laser light (tuned 
within realistic experimental conditions) is shown to induce an almost sixfold increase of the atom-mirror entan- 
glement and to be responsible for interesting dynamics between such mesoscopic systems. In order to assess the 
advantages offered by the proposed control technique, we compare the time-dependent dynamics of the system 
under constant pumping with the evolution due to the modulated laser light. 



Over the last few years we have witnessed a constant-pace 
advance in the control of mesoscopic systems at the quan- 
tum level. Examples range from quantum interference experi- 
ments with large organic molecules yj to the preparation and 
detection of multi-particle entanglement in superconducting 
devices All these efforts are contributing very signifi- 
cantly to the exploration of the quantum-to-classical boundary 
and the ascertainment of the existence of actual limitations in 
the enforcing of genuinely quantum mechanical behavior in 
mesoscopic systems. From a more applied viewpoint, these 
advances contribute to the quest for the manipulation of infor- 
mation encoded in open systems operating beyond the micro- 
scopic scale. 

Cavity quantum optomechanics has preponderantly 
emerged, recently, as a very interesting arena for the study of 
quantum features at the meso-scale. Its aim is to control the 
quantum state of mechanical oscillators by their coupling to 
a light field |3]. Recent advances in this context include the 
realization of the quantum-coherent coupling of a mechanical 
oscillator with an optical cavity 101 where the coupling rate 
exceeds both the optical and the mechanical decoherence rate, 
and the laser cooling of a nanomechanical oscillator to its 
ground state [jst]. When microwave radiation is used instead 
of optical fields, an analogous sideband cooling of the motion 
of a microscopic cantilever has been demonstrated, leading to 
a mechanical thermal occupation of ~^1 |@]. 

Such new technologies open unprecedented possibilities in 
the design of hybrid quantum architectures whose elemen- 
tary building blocks are physically implemented by systems 
of different nature. Achieving accurate control at their in- 
terfaces is quickly becoming of technologically and founda- 
tionally paramount importance. In this respect, the coherent 
coupling of a superconducting flux qubit to a spin ensem- 
ble has been recently reported [[3, while architectures for the 
interaction between ultracold atoms and mechanical systems 
have been demonstrated fsl. Furthermore, recent experiments 
have reached the regime of strong coupling between a Bose- 
Einstein condensate (BEC) and an optical cavity ||^[l3l, sim- 
ulating optomechanical effects where the mechanical system 
is embodied by quantum phononic waves of the BEC. 

Here we consider a hybrid situation where we combine 
technology coming from the traditional optomechanical set- 



ting with the potential of the BEC-cavity experiments and 
proposing a new BEC-optomechanics apparatus. The orig- 
inal idea ifTHl is founded on the insertion of a BEC in an 
externally driven optical cavity whose end-mirror oscillates 
around an equilibrium position. For a small number of me- 
chanical phonons, the composite system made of the cavity 
field, the BEC phononic mode, and the vibrating mirror is 
endowed with genuine multipaitite entanglement that can in 
principle be revealed by all-optical measurements lfl2l . We 
propose the active driving control of such a hybrid system re- 
alized by time-modulating the intensity of the driving field. 
Inspired by recent works on the control of optomechanical de- 
vices ifTsilTill . we show that by using a monochromatic mod- 
ulation, entanglement between two mesoscopic systems, the 
mirror and the BEC, can be created and controlled [15] . Fur- 
thermore, by borrowing ideas from the theory of optimal con- 
trol lfl6ll we show that, with respect to the unmodulated case, a 
sixfold improvement in the degree of generated entanglement 
is in order We interpret such performance in terms of the oc- 
currence of a special resonance at which the building up of 
entanglement is favored. Our results demonstrate the viabil- 
ity of the optimal control-empowered manipulation of open 
mesoscopic systems for the achievement of strong quantum 
effects, even in the hybrid context, of which the system that 
we study is a significant representative. 



I. THE PHYSICAL MODEL 

We start describing the hybrid optomechanical setup at 
hand lfTll[T2ll . which consists of a Fabry-Perot cavity with a 
vibrating end-mirror The cavity is pumped by a laser that 
couples to the mechanical mirror and an intracavity BEC. The 
Hamiltonian of the system (in a frame rotating at the fre- 
quency lol of the pump field) reads 



n = nc + nA + nM + Uag + Umc- 

The Hamiltonian of the mirror is 
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where m is the effective mechanical mass, is the free me- 
chanical frequency and q (p) is the position (momentum) op- 
erator of the mirror The Hamiltonian of the driven cavity is 

Hc" = fi(^c^^L)a^a~ift?7(a ~ a^), (3) 

where wc is the cavity frequency, a is the cavity field's anni- 
hilation operator, and 77 = y/2KTZ/liujL accounts for the laser 
pumping [TZ is the laser power and k is the cavity decay rate] . 
The BEC is taken to be weakly interacting, allowing the sepa- 
ration of the atomic field operator into a classical component 
(the condensate wave function) and a quantum one (the fluctu- 
ations), expressed in terms of Bogoliubov modes. The cavity 
is strongly coupled to the mode with wavelength Ac/ 2 where 
Ac is the cavity-mode wavelength |l^|Tl|]. We call uii, the fre- 
quency of the Bogoliubov mode and c (c^) the corresponding 
annihilation (creation) operator As the condensate is at low 
temperature, any thermal fluctuation of the atoms is negligi- 
ble. The free Hamiltonian of the Bogoliubov mode is given 
by Ha = hojbc'c, while the atom-cavity coupling is jgl fTTIl 



Hac 



where g is the atom-cavity coupling strength. No is the con- 
densate population, and is the atom-cavity detuning. In 
the second term, the atom-cavity coupling rate is denoted 
by The complete derivation of its form can be found in 
Ref. ifrHl . The position and momentum quadratures of the 
Bogoliubov mode are Q={c + c')/^/2 and P=i[& — c)/ 
respectively. Finally, the mirror-cavity interaction is given by 
T~Lmc^ — hxqa^a with x — ^c/L the mirror-cavity cou- 
pling rate (L is the length of the cavity). As can be seen 
by comparing T-Lmc with the second term of Eq. (|4]i, the 
BEC dynamics is analogous to a mechanical oscillator un- 
der the action of radiation pressure. As no direct coupUng 
term T-Lam is present in Eq. ([T]i, the interaction between the 
atomic mode and the mirror is mediated by the cavity. The 
relevant degrees of freedom of the system are grouped in the 
vector (f)^ = (x, y, q,p, Q, P), where the cavity position and 
momentum-like quadrature operators are x 



and y 



1) I \f2, respectively. Under intense laser 



pumping the operators can be linearised and expanded around 
their respective classical mean values 0^,1 such that 0; ^• 
(\)s^i + (50i, where (50-^ = (5y, 8q^ 6p, 6Q, SP) is the vec- 
tor of zero-mean quantum fluctuations for each operator in 
(f>. Here, the mirror position and momentum operators have 
been rescaled to dimensionless quantities as g = h/mumq, 
and p = \J hmujraP- In the hybrid optomechanical system 
considered here, two sources of noise should be taken into ac- 
count. The first comes from photons leaking out of the cavity, 
while the second is due to the mechanical Brownian motion 
performed by the mirror, which is typically in contact with a 
thermal bath at temperature T. The open-system nature of the 
problem at hand allows for the establishment of a stationary 
state. In fact, the classical values (^s,i can be calculated by 
solving the steady-state Langevin equations 117!], which leads 
to the new equilibrium positions for the mechanical mirror and 



the harmonic oscillator embodied by the Bogoliubov mode 
= ^^\"-f\ and, Qs = — -^-Ls^. Here we have introduced 

the mean intracavity field amplitude |asP = 77^/(A^ + k^) 
and the total cavity-pump detuning (modified by the radiation 
pressure mechanisms and the shift of the cavity frequency due 
to its off-resonant coupling with the atomic medium) 



A 



U}c - i^L + 



2Aa 



CQs 



(5) 



As the last three terms in A are typically very small compared 
to the bare detuning ujc—^^l, we neglect any bistability effect 
and assume A to be an independent control parameter 

As for the fluctuations, their dynamics can be described by 
the following vector equation 



(6) 



where we have introduced the input-noise vector AA^ = 
{\/2)i6xin, \/2fiSyin, 0, ^/y/hmujm, 0,0). The drift matrix 
(4) /C, given below, depends on the scaled coupling parame- 

and the mirror dissipation rate ■j=uj,n/Q 



ter x^x'\/Wrn^r, 
[where Q is the mechanical quality factor]. 
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The photon-leakage from the cavity is accounted for by the 
input-noise operators 



Sxi 



6 a, 



5 Hi, 
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and = i ^"™ 
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with {Sa,n) = (Sain) = and (5a™(i)M|„(i')> = Sit - t'). 
The Langevin force operator ^ in A/^ models the effects of 
the mechanical Brownian motion. In the limit of a high me- 
chanical quality factor, such noise can be faithfully consid- 
ered as Markovian, as entailed by the asymptotic form of 
the auto-correlation function {^{t)^{t'))2±h'ym/ /3BS{t—t'), 
where f3B=fi/2kBT and ks is the Boltzmann constant. 

We consider a viable detection scheme to observe entan- 
glement between the various bi-partitions. Our proposal is 
linked to the method put forward in Ref. iT2|| , i.e. on the map- 
ping of atom-field or mirror-field entanglement (and thus the 
atom-field one) into fully accessible all-optical quantum cor- 
relations by means of two extra fields that interact (locally) 
with the abovementioned subsystems. Both the operations are 
within reach of state-of-the-art experiments and, in fact, have 
been recently implemented ifisi \l9A with high efficiency. Our 
revelation scheme will thus be affected by the limitations of 
such methods, which embody the forefront of weakly disrup- 
tive detection schemes in such mesoscopic scenarios. 



3 




FIG. 1: (Color online) (a) Entanglement Em Ait) between mir- 
ror and atoms against Hit and A/tJm- (b) Same as panel (a) for 
EcM{t) = EcA{t). Parameters: bJm/2-K = 3 x 10^ s"^ 
T = 10/iA'; = 3 X 10"*; m = 50 ng; 7^ = 50 mW, cavity 
finesse F = 10*, and C = cavity length is L = 1 mm from 
which K = ■nc/2LF (c is the speed of light). All plotted units are 
dimensionless. 



II. ENTANGLEMENT DYNAMICS 

In Ref. |fl2ll . the stationary entanglement within the hybrid 
optomechanical system has been considered. Here we focus 
on the dynamical regime where the evolution of the entan- 
glement is resolved in time. We consider the fully symmet- 
ric regime encompassed by equal frequencies for the Bogoli- 
ubov and mechanical modes (i.e. ~ tUm) and identical 
coupling strengths in the bipartite cavity-mirror and cavity- 
atom subsystems (that is, we take C — x). We are partic- 
ularly interested in the emergence of atom-mirror entangle- 
ment at short interaction times. The analysis conducted in 
Ref. lfl2ll has shown it to be absent for t — > oo. The linear 
nature of Eq. (|6j preserves the Gaussian nature of any initial 
state of the overall system. This allows us to fully charac- 
terise the entanglement evolution through the covariance ma- 
trix = (5(/ij}) — (5(/)i)((5^j). Using this definition 
and Eq. (|6]l, the dynamical equation regulating the evolution 
of the covariance matrix can be written as 

V = /CV + V/C^ + 2?, (9) 



where we have introduced the noise matrix V = 
diag[K, K, 0, 7(2n+l), 0, 0] with n = [exp(2/3BW,„) — 1]^-^ 
as the thermal mean occupation number of the mechanical 
mode. In Eq. ^ we assume that the mean values of the me- 
chanical and optical quadratures reach their stationary values 
much faster than the fluctuation dynamics (this is always ver- 
ified in our calculations). Such an inhomogeneous first-order 
differential equation is solved assuming the initial conditions 
V{Q) = diag[l, 1, 2n+l, 2n+l, 1, l]/2, which describe the 
vacuum state of both the cavity field and the BEC mode and 
the thermal state of the mechanical system. Physicality of 
the covariance matrix has been thoroughly checked by consid- 
ering the fulfillment of the Heisenberg-Robinson uncertainty 
principle and checking that the minimum symplectic eigen- 
value V = mill cig(iwV) is such that > ^. Here used 
the 6x6 symplectic matrix uj = (Bj^ii<Jy with dy the y- 
Pauli matrix. The entanglement measure that we use here to 
quantify entanglement between any two modes a and /3 is the 
logarithmic negativity 1I20I1 , defined as E^p — — log|2;^mi„|, 
where i^mm is the smallest symplectic eigenvalue of the ma- 
trix Vj^ = PVapP, for a,l3 = C,A,M [the latter being 
the labels for the cavity, atomic and mirror modes, respec- 
tively]. The reduced covariance matrix Vafs contains the en- 
tries of V associated with modes a and /? while, by inverting 
the momentum quadrature of /?, matrix P = diag(l, 1,1,-1) 
performs the partial transposition in phase space. The atom- 
mirror entanglement EMA{t), whose time evolution is shown 
in Fig. [lla) for different values of the effective detuning 
A, gradually develops and reaches its peak value as the 
cavity-atom and cavity-mirror entanglement drop to a quasi- 
stationary value. As no direct atom-mirror interaction exists 
in this system, mediation through the cavity mode essentially 
results in a delay: quantum correlations between the atoms 
(the mirror) and the cavity mode must build up before any 
atom-mirror correlation can appear This is clearly shown in 
Fig. Eb), where EcAi{t) — EcA{t) reach their maximum 
well before Em a (t) starts to grow. The atom-mirror entan- 
glement is non-zero only within a very short time window, 
signalling the fragility of quantum correlations resulting from 
only a second-order interaction between the BEC and the cav- 
ity end-mirror These results go far beyond the limitations of 
the steady-state analysis conducted in lfl2ll and prove the exis- 
tence of a regime where all the various reductions obtained by 
tracing out one of the modes from the overall system are in- 
separable, a situation that, within the range of parameters con- 
sidered in our investigation, is typical only of a time-resolved 
picture. 



A. Optimal control of the early-time entanglement 

We now consider the effects of time-modulating the exter- 
nal pump power TZ, which is now considered a function of 
time. In turn, this implies that we now take 1] ri{t) and 
study the time behavior of the entanglement Em a set between 
the atoms and the mirror We will show that a properly opti- 
mised 7]{t) can increase the maximum value of EMA{t) for 
values of t within the same time interval t where atom-mirror 
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entanglement has been shown to emerge in the unmodulated 
case. We assume to vary ri{t) slowly in time, so that the clas- 
sical mean values (f>s adiabatically follow the change in rj{t). 
This approximation is valid as long as the number of intra- 
cavity photons is large enough to retain the validity of the 
linearization procedure and the time-variation of ri{t) is slow 
compared to the time taken by the mean values to reach their 
stationary values. For all cases considered here we have ver- 
ified the validity of such assumptions. The dynamics of the 
covariance matrix is thus still governed by Eq. (|9) with the 
replacement /C — > /C(t). In the following, we use the value 
of A = 2.7u!,n which maximises the short time entanglement 
Ema- 

Inspired by the techniques for dynamical optimization pro- 
posed in lfl6ll . we call 770 the unmodulated value of i] and take 



r/{t) = ryo + ^ [^j cos{ujjt) + Bj sin(a;jt)] , 
i=i 



(10) 



where coj ~ 27rj/r + Sj are the harmonics and Sj is a small 
random shift. The coefficients Aj and Bj are chosen in a way 
that the total energy brought about by the time-modulated field 
is the same as the one associated with the unmodulated case. 
We then set the time-window so that t = 3.4k^^, when we 
observe the maximum value of Ema in the unmodulated in- 
stance. The other parameters are as in Fig. [T] We then look for 
the parameters Aj and Bj that maximise the value of Ema (t) 
for a given set of random shifts 6j . We use standard optimi- 
sation routines to find a (local) maximum of Ema{t)- We 
repeat the search of the optimal coefficients for different val- 
ues of 5j and take the overall maximum. The corresponding 
results are presented in Fig. |2] where we show the optimal 
modulation 77 (t) and the optimal EMA{t). These findings are 
also compared to the case without modulation. The maxi- 
mum value attained in the interval [0; r] is EMA{t) — 0.05 
which is 2.5 times larger than the case without modulation, 
thus demonstrating the effectiveness of our approach. 



Fig. |3] (a) (inset) revealing a sharp resonance with a maxi- 
mum value of Ema — 0.12 for S = S ^ 0.79k (no fur- 
ther peak appears beyond this interval). This arises as a result 
of the effective interaction between atoms and mirror medi- 
ated by the cavity field. As shown in Fig. [TJa) the entangle- 
ment dynamics strongly depend on the effective detuning giv- 
ing rise to such optimal behavior Similar results have also 
been observed in Ref [| 1 3il . The analysis of the evolution of 
EcM{CA) and Ema, shown in the main panel of Fig. [3] (a), 
reveals that while £'ca/(ca) develops very quickly due to the 
direct cavity-atoms and cavity-mirror couplings, Ema grows 
in a longer time lapse, during which the cavity disentangles 
from the dynamics. The quasi-asymptotic value achieved by 
Ef^iA reveals a sixfold increase with respect to the unmodu- 
lated case. This behavior is worth commenting as it strength- 
ens our intuition that any atom-mirror entanglement has to re- 
sult from a process that effectively couples such subsystems, 
bypassing any mechanism giving rise to multipartite entangle- 
ment within the overall system. Finally, we discuss the results 
achieved in the long time case by adopting an optimal-control 
technique similar to the one used for enhancing the short time 
entanglement. We have considered the periodic modulation of 
the intensity at the frequency E given by 



ri{t) 



1 - ^ A„ sm{nT,t) + B„ cos{nT,t) 



(11) 

and looked for the coefficients {An, -B„} optimizing EMA{t) 
at long times with the constraint: J2n=T i^n + ^n) — 1' 
ensuring that no instability is introduced in the dynamics 
of the overall system. In our simulation, rimax = 8 has 
been taken to limit the complexity of the modulated signal. 
The resulting optimal coefficients give the periodic modu- 
lation shown in Fig. [3] (b) and a maximum entanglement 
Ej^jA — 0.17 which is about 30% larger than the results ob- 
tained for the monochromatic modulation. This demonstrates 



B. Periodic modulation: long time entanglement 

In the two situations analyzed so far [constant laser inten- 
sity 7] and optimally modulated ri{t)], the atom-mirror entan- 
glement Ei\iA{t) is destined to disappear at long times. A 
complementary approach based on a periodic modulation of 
the laser intensity ri{t) was used in the pure optomechan- 
ical setting lfl3ll to increase the long-time light-mirror en- 
tanglement. Here we use a similar approach by assuming 
the monochromatic modulation of the laser-cavity coupling 
ri{t) — t/q + rjQ [1 — siii(I]t)], where S is the frequency of the 
harmonic modulation, 77Q — Arj^ = rio/2, 770 being the same 
constant coupling parameter taken before. These choices en- 
sure that the approximations used in the dynamical analysis 
are valid. After the transient dynamics, the covariance matrix 
and, in turn, Ema (t) become periodic functions of time. In 
order to achieve the best possible performance at long times, 
we compute the maximum of Ema (t) after the transient be- 
havior, scanning the values of S. The result is shown in 




FIG. 2: (Color online) Entanglement dynamics Em Ait) (dashed 
line) with the optimal laser intensity modulation (solid line). 
The entanglement EMA{t) for constant rj{t) = r;o is also shown 
(dashed-dotted line). All plotted units are dimensionless. 
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FIG. 3: (Color online) (a) Dynamics of the cavity-mirror and cavity- 
atoms entanglement Ecm,ca (solid line) and atoms-mirror entan- 
glement Em A (dashed line) for E = E ~ 0.79k. Inset: Maximum 
Em A for long times with a periodic modulation as a function of the 
frequency E. (b) Optimal periodic modulation ri{t) (solid line) for 
one period of time 27r/E compared to the monochromatic modula- 
tion (dashed line). All plotted units are dimensionless. 

the powerful nature of our scheme. Our extensive muhiple- 
harmonic analysis is able to outperform quite significantly 
the single-frequency driving scheme addressed above and dis- 
cussed in Ref. |[l3|, proving strikingly the sub-optimality of 
the monochromatic modulation, both at short and, surpris- 
ingly, at long times of the system dynamics. This legiti- 
mates experimental efforts directed towards the use of time- 
modulated driving signals for the optimal control of the hybrid 
mesoscopic systems addressed here and similar ones based, 
for instance, on the use of a vibrating membrane ll2l1l or a 
levitated nano-sphere ll22ll instead of the BEC. 



It should be noted that an adiabatic approach is used in 
Ref. ifisll to find an effective Hamiltonian. When performed 
in our hybrid optomechnical scheme, such technique would 
remove the dynamics of the Bogoliubov modes of the atomic 
subsystem. 

We consider the robustness of our protocol with respect 
to inaccuracy in the control of the value x = C follows: 
after finding the best periodic modulation assuming x = C 
we ran again the simulations with the same modulation with 
X — 1-1^ and X = 0.9^. These values correspond to a 10% 
inaccuracy in the nominal values of x and ^. We found that 
the maximum entanglement is at most only 3% less than the 
original value thus confirming the stability of our result. No- 
tice also that if the imbalance between x and ^ is known, for 
example by a calibration measurement, we can in principle 
run the optimization including the actual values of x and ^ 
therefore aiming at a larger entanglement value. 

III. CONCLUSIONS 

We have demonstrated that the modulation-assisted driving 
of a hybrid optomechanical device gives rise to interesting 
and rich entanglement dynamics, surpassing the limitations 
associated with a steady-state analysis and a constant pump. 
A significant improvement in the genuinely mesoscopic en- 
tanglement between the mode embodied by the atomic sys- 
tem and the mechanical one can be achieved by pumping the 
cavity with a modulated driving field, in both the short-time 
case and the long-time case. We have shown the existence of 
modulations that are able to beat the performance of simple 
monochromatic driving in terms of the maximum entangle- 
ment created between the mirror and the atomic system, thus 
favoring the creation of entanglement at long times. 

Our study strengthens the idea that important advantages 
are in order when optimal control techniques are implemented 
in the open-system dynamics of mesoscopic devices. This 
contributes to the current quest for the grounding of such 
approaches as valuable instruments for the control of meso- 
scopic (and multipartite) systems of various realizations, in- 
cluding interesting configurations of current experimental in- 
terest 1I2TII22I1 to which our framework can be fully applied. 
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